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Given the state of a quantum system, one can calculate the expectation value of any observable 
of the system. However, the inverse problem of determining the state by performing different 
measurements is not a trivial task. In various experimental setups it is reasonably straightforward 
to reconstruct the state of a quantum system employing linear tomographic technique. In this way 
the elements of the density matrix can be linearly related to a set of measured quantities. But 
since different observables of a quantum system are not commuting with each other, one often has 
to perform series of successive measurements of observables which cannot be done simultaneously. 
Simultaneous measurement of observables cost less time and energy and is more beneficial. In 
this paper we review the strategy of quantum state tomography with simultaneous measurement of 
commuting observables. This can be done by introducing an assistant system of which the state 
is known. We show that the interaction between the assistant and the system of interest within 
different frame works allows the reconstruction of the state of the system. Specifically, we consider 
a two-level system and reconstruct its initial state by introducing an assistant which can be either 
another two-level system or a single cavity mode of the electromagnetic field. 

Keywords: quantum state tomography, simultaneous measurement of non-commuting observables, Jaynes- 
Cummings model, state determination, coheren state of light, two-level system, spin-Tj system 
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I. INTRODUCTION 

In classical physics, the state of a system is character- 
ized by specifying the values of all physical quantities, 
for instance the positions and the velocities of the parti- 
cles that constitute the system. In quantum mechanics 
the situation is complicated by the fact that the physi- 
cal quantities are mathematically represented by specific 
type of operators called observables, which in general are 
elements of a non- commutative algebra. Hence their val- 
ues cannot be simultaneously specified, as emphasized in 
the Heisenberg's uncertainty principle. Instead, the mea- 
surement results of each observable is characterized by a 
probability distribution, which involves statistical fluctu- 
ations. The "state of the system" is then represented by 
a mathematical notion that allows us to express the prob- 
ability distribution of all the observables for an ensemble 
of identically prepared systems. 

This is best described in the statistical interpretation 
of quantum mechanics^—, to which we adhere. In this 
interpretation, the state of a system is a mathematical 
object from which we can derive any probabilistic predic- 
tion about the physical quantities attached to this sys- 
tem. One typically imagines some experimental appara- 
tus and procedure which "prepares" this quantum state; 
the mathematical object then reflects the setup of the 
apparatus. This way, the quantum state accounts for the 
full information available about the preparation of the 
system, from which we wish to derive consequences for 
future experiments. Since this knowledge is probabilis- 
tic it does not refer to a single system or a single event. 
Thus what we call a state, which is most of the time a 
mixed one, characterizes a statistical ensemble of systems 
of the same type, which are all prepared under identical 
physical considerations. The state is thus a mathemati- 
cal representation of the result of a certain state prepa- 



ration procedure; it accounts for our information about 
this preparation and upon knowing it we can elaborate 
consistent probabilistic predictions. A standard tool to 
implement the statistical definition of state is the den- 
sity matrix, which generalizes the pure state represented 
by a wave function. Indeed, there is no conceptual dif- 
ference between wave function and density matrix which 
are both mathematical means for evaluating expectation 
values of the observables of the system or probabilities. 

In the frame work of the statistical interpretation, the 
laws of quantum mechanics can be summarized as fol- 
lows^— : 

• An observable O is represented by a self-adjoint 
linear operator acting on the Hilbert space pertain- 
ing to the system. It has a spectral representation, 
O = ^iOiPi where Oi are the eigenvalues of O 
and Pi are the orthogonal projection operators re- 
lated to the orthonormal eigenvectors of O, i.e., 
Pi = ^2 m \m, Oi)(m, Oi\. The parameter m labels 
the degenerate eigenvectors of O. 

• The state of a system at a given time is represented 
by its density matrix, p, which is a self-adjoint op- 
erator in Hilbert space with a unit trace. The den- 
sity matrix should also be semipositive to ensure 
that any variance of the observables of the system 
is non-negative. Pure states correspond to the spe- 
cial case 



p = p. 



(1) 



• The dynamics of the system can be obtained by 

p{t) = u(tM)p{h)u\tM), (2) 

where [/(£, to) is the unitary time evolution opera- 
tor. 



• Given the density matrix p of a system, one can 
find the expectation value of any observable O of 
the system in the considered situation as 

(6) = tT[p6), (3) 

where tr[- • • ] stands for the trace of a matrix. 

Let us emphasize that through out this paper the opera- 
tors are always distinguished by a " sign. 
Now we wish to face an inverse problem. Consider an en- 
semble of systems S prepared in some well-defined, but 
unknown, fashion. Nothing is a priori known about their 
state p, in other words, the probability to observe some 
result or another in the measurement of an observable is 
unknown. The following question then is of our interest. 
How can one determine the density matrix of this ensem- 
ble by identification of a set of commuting observables, 
the measurement of which permits the precise determi- 
nation of pi In other words, how can one determine the 
quantum statistical operator that describes the prepara- 
tion of the system? 

Procedures of reconstructing the quantum state from 
measurements are known as quantum state tomography. 
Recently, they have found some applications in quantum 
information processing^. For example, in quantum cryp- 
tography one needs a complete specification of the qubit 
state both as it is emitted from the source and as it is 
received after transmission^. 

In the simplest example of a spin-^ system or equiv- 
alently any two-level quantum system the state is de- 
scribed by a 2 x 2 matrix. In the two-dimensional Hilbert 
space, any observable is a linear combination of Pauli op- 
erators, which satisfy 



a a = 1, a = x,y,z, 
and arc represented by the Pauli matrices 



(4) 



1 

1 



-i 

i 



1 
-1 



(5) 



A state is characterized by three real numbers: one for 
the diagonal elements of the 2x2 density matrix p, and 
two for its off-diagonal elements. Equivalently, we can 
introduce the polarization vector, r, the components of 
which are the expectation values of the Pauli matrices. 



r a = tr [pa a ] , 



x,y,z. 



(6) 



Once we know the value of these parameters, we are 
able to determine the value of the density matrix, making 
use of the identity 



r ■ a 



(7) 



Thus, according to the above argument, one has to 
perform three incompatible measurements for the un- 
known state determination, e.g., measuring the spin com- 
ponents along the x-, y- and z- axes via Stern-Gerlach 
setup. However, during the measurement procedure of 
each component one looses the information about the 
two other components, since the spin operators in differ- 
ent directions do not commute. Thus, to determine the 
state of a spin-i system, one needs to use three sets of 
Stern-Gerlach measurements performed along orthogo- 
nal directions. In this approach, the state of any two- 
level system, represented by a 2 x 2 density matrix p, can 
be fully determined only through measurement of three 
linearly independent observables which do not commute 
and cannot be simultaneously measured. 

It has been proven that the unknown density matrix of 
such a system S, in particular the full polarization vector 
of a spin-i system, can be determined indirectly. This 
can be done by means of a single set of measurements 
performed simultaneously on S and an auxiliary system 
(assistant or ancilla). The assistant (or ancilla) starts its 
evolution from a known stated— 

The suggested strategy is the following: Initially S lies 
in the unknown state that we wish to determine, while 
the state of the assistant A is known. During some time 
lapse S and A interact in a known fashion. As a re- 
sult their joint state is modified: it involves correlations 
and keeps memory of the initial state of S. Two com- 
muting observables of the combined system (system + 
assistant) are then simultaneously measured. Repeating 
this process provides then the necessary statistical data: 
the expectation values of the observables and also their 
correlation. We will show that one can infer the three 
components of the initial polarization vector of S, and 
hence the state of the system from these three sets of 
data. This type of information transfer is remarkable 
because initially an unknown information was embedded 
in the matrix elements of p, or equivalently in the com- 
ponents of the polarization vector of S. It had a quan- 
tum nature, and could not be represented by an ordinary 
probability distribution, due to the non-commutation of 
the three components of the spin operator. After the 
interaction between S and A this unknown initial infor- 
mation about S, together with the known information 
about A, is redistributed among the matrix elements of 
the joint density matrix of the overall system, S+A. How- 
ever, the resulting classical joint probability distribution 
for the observables of the system S and the assistant A 
can keep full memory of the initial quantum informa- 
tion about S. The process on which we rely, amounts 
to a transformation of quantum information into classi- 
cal information, which can be gained by a classical type 
of measurement involving commuting observables only. 
This measurement modifies the state of S+A, but it can 
preserve all the matrix elements of p. The idea of trans- 
forming quantum into classical information by using an 
assistant system A was first proposed by D' Ariano^ who 
showed the possibility of mapping the density matrix of S 



onto a single observable of S+A. It was explicitly imple- 
mented in a dynamical form by Allahvcrdyan et al. — . In 
particular, the authors^ showed that one can determine 
the unknown state of a spin-i system with a single ap- 
paratus by using another spin-i assistant. This idea was 
recently implemented by Peng et al£ who used pulses to 
induce the proper dynamics of the interaction between 
the spin- 2 system and its assistant. They verified the 
initial state of the system obtained from this procedure 
with the result of direct measurements of the three com- 
ponents of the spin vector of the system. 
Aquino et al*& and Mchmani et alM- showed that the ini- 
tial state of a two-level system can be characterized by 
introducing a single mode of a coherent light instead of 
introducing another spin, so the dimension of the Hilbert 
space of the assistant can be larger than that of the sys- 
tem. More specifically, Aquino et al. showed that in 
the resonant case, i.e., when the frequency of the two- 
level system is equal to the cavity mode frequency, it 
is not possible to determine the initial state of a two- 
level system by measuring the the energies of the system 
and the assistant. However, this scheme is still applica- 
ble for spin- 2 systems because in order to recover the 
initial state of the spin-i system one can measure the 
transversal component of spin [x or y component) and 
the number of the photons in the cavity^. In general, 
for two-level systems other than spins, it is rather diffi- 
cult to measure the transversal component since it cannot 
be defined well. While the z-component of quasi-spin is 
related to the level occupation and thus the energy of 
the two-level system. The crucial point in this case is 
that there should be a detuning between the frequency 
of the field and that of the system of interest in order 
to invert the relavent relations between simultaneously 
measured observables in one hand and the elements of 
the initial density matrix of the two-level system on the 
other hand&. 

This overview is organized as the following: In sec- 
tion |TT] we review the idea of determining the state of a 
spin- 2 system by letting it interact with another spin- -^ 
system which starts its evolution from either a completely 
disordered state or a pure state. Section IIIII is devoted 
to the situation where the assistant is considered to be 
a coherent single mode of the electromagnetic field. The 
last section IPvH presents our conclusion. 



II. SPIN AS AN ASSISTANT 

As it was mentioned before, the aim is to find an in- 
direct procedure to determine the unknown state of a 
quantum system. The desired procedure only consists 
of measurements of commuting observables, which there- 
fore can be performed by means of a single apparatus. 
For the quantum system we consider a two-level system, 
S, the state of which we wish to determine. Let the sys- 
tem S be coupled to an auxiliary two-level system A. 
A is in a known state. At later time we measure the 



value of one observable of S and one of A. Since the two 
observables belong to two different Hilbert spaces, they 
commute and therefor the measurement process can be 
performed simultaneously. Moreover, the correlator of 
two observables can be read off from coincidences. 

Let the initial unknown state of S be a density matrix 
of the form 



P 



\{ 



1 



where 



def 



tr [p &i 



x,y,z. 



(8) 



(9) 



The state is called pure if \r\ = 1, in which case the 
eigenvalues of p are and 1. \r\ < 1 represents a mixed 
state, and \r\ > 1 is physically excluded. 

We choose the state of the assistant, represented by R 
as 



R = - (1 + \s z 



0< A< 1 



(10) 



where s x , s y , and s z arc the Pauli matrices in the Hilbert 
space belonging to the assistant A. 

Initially there is no interaction between S and A. There- 
fore the initial state of the overall system, Oo, can be 
written as 



n = Rp 



(l + A)(l + f-<7) 

(l-A)(l+f-<7) 



(11) 
Now we let the two systems interact for some time. The 
interaction can be described with the help of a 4 x 4 



unitary matrix U 



-iH 



where we set t = 1. Here, 



we don't specify our Hamiltonian and consider a general 
unitary matrix and we parametrize it such that it gener- 
ates a proper time-evolved overall density matrix at later 
time t = 1, given by Clf such that the initial state of S 
can be read off easily. The observables of which the mea- 
surements yields the determination of the initial state of 
S are the final polarization of each spin of the overall sys- 
tem S+A£. They can be measured simultaneously and 
the correlator of the two can be derived from the gath- 
ered data. We show that it is possible to read off the 
initially unknown state of the system S from the three 
above mentioned sets of data. Let us decompose U into 
the following 2x2 block matrix, 



U = 



A C 

b b 



(12) 



and express the unitarity of U in terms of the 2x2 ma- 
trices A, B, C, D in the Hilbert space of S. The polar 
decomposition of A and B yields 



A = vk. 



B = wk\ 



(13) 



where v, and w are unitary 2x2 matrices while k and k' 
are semi-positive Hcrmitian 2x2 matrices. Since v, and 



w are unitary, it is easy to see that k and k are the non- 
negative square roots of A* A and W B, respectively. If k 
and k! have a vanishing eigenvalue, these representations 
of A and B still hold but are no longer unique. We shall 
restrict ourselves to the case where k and k' are strictly 
positive. 
The condition UW = 1 implies 

C& = 1- AA\ f)b ] = 1-BB\ (14) 
iS f + CD* = 0, (15) 

while tnU = 1 implies 

#4 + ^5 = 1, C^ + b ] D = 1, (16) 
i t C' + i? t i) = 0. (17) 

Implementing (|16p on the polar decomposition of A and 
5 given by (|T5)) yields 



fc' = \/l-fc 2 . (18) 

Thus CCT and Z5.D' can be simplified as 

C& = vk' 2 v\ DD* = wk 2 wK (19) 

Since k and k' are strictly positive and u and w are uni- 
tary matrices, we can define unitary matrices x and y 
such that C and D have the form 



C = iik'x, 



D = wky. 



(20) 



The remaining unitary condition A'C + B'D = reads 
kk (x + y) = 0. Again, since k and k are strictly positive 
this implies y = — x, which fixes y in a unique way. 
The unitary matrix U then becomes: 



U 



v 



w \k' 



k k' 



1 
x 



(21) 



In order to get a more symmetric form for U, we intro- 
duce a unitary matrix X such that (X*) 2 = x and define 
the matrices V = vX* and W = wX\ K = XkX\ 
K' = Xk'XK We can then write the 4x4 unitary trans- 
formation operator U as 



U 



V 

w 



K K' 

k' -k 



X 

o i"f 



(22) 



in terms of the three unitary matrices V, W, X and the 

non-negative hermitian matrices k and K' = \/l — K 2 . 
Since K and K' arc strictly positive, this decomposition 
is unique, provided we fix the signs of the eigenvalues of 
X = vff by some convention, for instance, (X+X*) > 0. 
Having the unitary matrix given by (|22[) we can calcu- 
late the state of the overall system at later time as 



flf = Utt U\ 



(23) 



where Ctg is given by ([TTj) . The obscrvables that can be 
simultaneously measured by means of the same appara- 
tus on the state £7/ are the z-componcnts of each spin. 
This corresponds to the following averages 



&) = 


= tr 


ClfS z 


) 


(^) = 


= tr 


n,fa z 


* 


Sz&z) '- 


= tr 


% (s z a z ) 



(24) 



We notice that the correlator (flf (s z a z )) can be 
recoverd form the gathered data of a z and s z via the 
number of coincidences. 

Inserting the unitary time-evolution operator in the 
expectation values (f2~4"| we get a linear relation between 
the gathered data of measurements of s z , a z and their 
correlator on one hand and the elements of the initial 
density matrix of S given by r x , r y and r z on the other 
hand. But before calculating the above mentioned 
expectation values, let us first parametrize the unitary 
time evolution operator U. 

Since K is a Hermitian matrix with < K < 1, we 
can parametrize it as 



K = cos 6 cos <f> + sin 6 sin 4>(x- a), 



(25) 



where % is a unit vector and 0<<f><8<^ — <p. It is 
straightforward to see that K' is given by 



K' = sin 9 cos (f> — cos 9 sin <fi(x ■ a), 



(26) 



Since the initial overall density matrix CIq is block di- 
agonal, multiplication of the unitary matrix X by a phase 
factor does not affect UfloW although it modifies U . 
Therefore, we can parametrize X as 



X = e^Cf*) = C o S Vj + j(|\ g) sin^, 



(27) 



where £ is a unit vector which we assume to be perpen- 
dicular to x for simplicity, and < ip < n. 
Parametrization of V and W can be done due to the fact 
that we are not interested in the off-diagonal block ele- 
ments of 51/. In other words, the three expectation values 
([Ml) do not require the determination of the off-diagonal 
elements of the overall density matrix and it would be 
sufficient to determine the action of V and W on the a z : 



y t CT z V r = ?7-CT, W**cr z W = ( 



a. 



(28) 



where r\ and C, are three dimensional unit vectors. 

Inserting the expression for flf from (j23|) into ([24]) us- 
ing the parametrization introduced by ([2"5 ]) -([2"8 ]l yields 

(s z ) = A cos 29 cos 2</> + X(x ■ r) sin 29 sin 20 cos 2ip 
+ [(l*x X) ■ H sin 26* sin 20 sin 2-0, (29) 



and 



((Tz) + (s z Vz) = A(% • ff) sin20 sin20 

+ (X ■ v)(x ■ r)(l - Acos20)(l - cos20)cos20 
+(x-?T) {£><X)-r (A-cos20)(l-cos20)sin20 

+(£*• ?7)(£ ' r)(cos20 + Acos20)(l - cos 20) 
+ (ff ■ r*) (cos 20 + A cos 20) cos 20 



(£ X ff) ■ r 



(A cos 20 + cos 20) sin 2-0. 



(30) 



Finally, (er 2 ) — {s z cr z ) can be obtained by transforming 
29 to 29 + ir and replacing ff with £ in (|3"0|) . 

For the sake of simplicity, we assume that £ is the unit 
vector in the a;— direction and that the unit vector x lies 
in the y— direction: 

(= (1,0,0), x= (0,1,0), fxx= (0,0,1). (31) 

Therefore the components of the two vectors ff and £ on 
£ and x can b e defined as 



def i» _ 


tie/ _ _ 

% = x-v, 


% = K x xl -^ 


c^r-c, 


C='x-C, 


C, = Kxx]-C- 



(32) 

Within the above choice of the unit vectors we can 
relate the measured values of the population difference 
of the two energy-levels of A and S to the initial state of 
the system S as 




= C 




(33) 



where C is a 3 x 3 matrix of coefficients the elements the 
elements of which are given by 



en 

Cl2 
Cl3 
C21 

<"22 



0, 



(34) 



= A sin 20 sin 20 cos 20, 

= sin 28 sin 2(f) sin 20, 

= {Vx + CO cos 20 + \{-q x ~C)cos20, 

= (7? y + Cy) cos 20 + A(?7 y — CO cos 20 cos 20 cos 20 

- sin 20 [A(7y z + C) cos 20 + (^ - G ) cos 20] , 
= ^{Vy + (y) sm 20 + (?7 y — GO cos 20 cos 20 sin 20 
+ cos 20 [(r) z + Q) cos 20 + \{r) z - G) cos 20] , 
= (Vx -G) cos 20 + \(t] x +( x )cos29, 

- (Vy ~ Cy) COS 20 + \{r} y + (y) COS 20 COS 20 COS 20 

- sin 20 [X(rj z - (z) cos 20 + {r) z + ( z ) cos 20] , 
= X(r] y — C y ) sin 20 + (rj y + Cy) cos 20 cos 20 sin 20 
+ cos 20 [(r} z - G) cos 20 + A(?y 2 + Q) cos 20] , 

(35) 

and T is the vector of constants given by: 



<"2:5 

C31 
C32 

C33 



cos 20 cos 20 
? A ( (r?j / -C J/ )sin20sin20 

Jj) y + Cy) sm 20 sin 20i 



(36) 



The elements of the initially unknown density matrix 
of S which are encoded by r are related to these expec- 
tation values, so they can be recovered if and only if the 
determinant of the coefficient matrix C is non-zero. With 
some algebra we can calculate the determinant of the co- 
efficient matrix, represented by D, as 

(1- A 2 ) sin 40 .. „, , „„, ^ 

— [(cos 20 + A cos 20) T) X C V 



8D 



sin 20 sin 20 2 

— (cos 20 — A cos 20) r} y C x ] 

+ ?7zCz;(cos20 — A cos 20) [A cos 20 

+ cos 20(A 2 cos 2 20 + sin 2 20)] 

— ^Cz (cos 20 + A cos 20) [A cos 20 

— cos20(A 2 cos 2 20 + sin 2 20)]. (37) 

Thus the initial state of the system S can be deter- 
mined from (s z ), (a z ) and (s z a z ) provided that the de- 
terminant D is non-zero. 

In what follows we consider two limiting cases: i) when 
the assistant is initially in a completely random state 
(A = 0), and ii) when it starts its evolution from a pure 
state, i.e. A = 1. Then we maximize the value of A over 
the parameters of U and reconstruct the initial state of 
S. We also try to understand what kind of spin-spin in- 
teraction Hamiltonian will yield the maximum value of 
D in each case. 



A. Assistant in completely disordered initial state 

Inserting A = in the expression for the determinant 
given by (|3"T|) yields 

D = — sin 20 sin 40 sin 20 [cos 20 (^ C y — r] y Cx) 

+ cos20sin20feCx + ??xG)]- (38) 

It is clear that this determinant is maximized over the pa- 
rameter if = ± -| . Furthermore, the maximum of D = 
jg sin 20 sin 20 [cos 20 (j] x C y — T] v Cx) + cos 20 sin20(r/ z G + 
VxCz)] over ff and C is reached when 



*-?-&°-ji>- 



Thus we have 
1 



D = — sin 20 sin 20\/ 1 - sin 2 20 sin 2 20. 
16 



(39) 



(40) 



The determinant (|40[) reaches its maximum value 1/32 
for = 7r/8 and = 7r/4. Such non-zero determinant 
guarantees the procedure of inversion and characterizing 
the initial state of the system. Inserting the above values 
of the parameters in the expressions for the expectation 
values of the z-component of spins, ([2T))) . (|30|) we can 
reconstruct the initial density matrix of S: 



r x = 

r y = 

r z = 2(s 



2(cr z >, 

-2(s z a z 



(41) 



We see that for a suitable choice of the evolution oper- 
ator U it is possible to determine the initial state of a 
spin-1/2 system implying an assistant which is initially 
in completely disordered state. 



1. Construction of a Hamiltonian 

In order to build up a Hamiltonian such that the value 
of the determinant \D\ is maximized, we consider the 
following feasible criteria: 

• The interaction between A and S is described by 
an anisotropic Hcisenbcrg Hamiltonian 



int — E ^*^**»j 



(42) 



where J Xl J y and J z are the couplings for x, y and z 
components. This type of Hamiltonian is employed 
to describe interacting spin systems in different ar- 
eas such as spin chains, nuclear magnetic resonance 
(NMR), quantum optics, spintronics, etc. 

• The possibly external fields acts symmetrically on 
A and S. Thus the total Hamiltonian can be writ- 
ten as 

H = J2 W°i + Si) + H iat , (43) 

where hi are external fields in x, y and z direction. 

Now let us try the following form: 

H = h x (a x + s x ) + J x (cr x s x + a z s z ) + J y a y s y , (44) 

which involves a partially anisotropic Hcisenbcrg Hamil- 
tonian and an external field acting on the x— direction. 
This Hamiltonian can be directly diagonalizcd. Its eigen- 
values Ek are 

Ei = —2J X — J y , Ei = Jy, 

E 3 = J X - P, 

E 4 = J X + p, (45) 

where we defined 



P=yJ(J x -J y ) 2 +4hl. 



(46) 



The corresponding orthonormal cigenstates \E k ) , k 
1, . . .4, read 

\Ei) = ^(0,-1,1,0), 

\E 2 ) = -±=(-1,0,0,1), 



\Eo 



1 



= (1,7,7,1), 



V2(l + 7 2 ) 

\E 4 )= -==L= =(l ;At ,/x,l), (47) 

V 2 (l + M 2 ) 



where 7 and fi are defined as 

7 

Now for t — 1 we can write 



J y — Ei 
2h x ' 

Jy — E 3 

2K 



U 



iH V^ 



-iE k 



E k )(E k \ 



(48) 
(49) 

(50) 



fc=i 



The expectation value of any observable O can be calcu- 
lated as 



(6) = tr 



tlntfldft 



E £ 

fc,J=l a=x,y,z 



j(E k -Ei) 



{E k \6\Ei){Ei\a a \E k )r a 



(51) 




where we have used (fTTj) . Calculation of (s z ), (& z ) and 
theire correlator yields 



(52) 



where V is a 3 x 3 matrix whose determinant is 

^ = (J x -Jyfhls^AJ x )s^^ (g3) 

The determinant (|53p reaches its maximum value of 1/32 
for 



■■' :r 



J ,, 



tt(1 - V8) 



± 



2-y/8' 



(54) 



Thus the optimal measurements for the completely ran- 
dom initial state of the assistant are reachable with par- 
tially anisotropic Heisenbcrg Hamiltonian and one ex- 
ternal magnetic field acting on both systems along the 
x— direction. 

We note that the physical reason of maximizing the de- 
terminant D is to suppress the effect of the errors made 
during the measurement process. 



B. Assistant in a pure initial state 

Considering the assistant starts its evolution from a 
pure state is equivalent to set A = 1 in the general ex- 
pression for the determinant given by (|37[) . which yields 

D = -sin26»sin2</>(cos 2 0- cos 2 6) (rj z ( x - rj x ( z ). (55) 



The maximum value of the determinant, \D\ = l/12-\/3, 
in this case is reached when fj and C, are perpendicular to 
each other 

4> = ±j, and sin 2 (29) = 1/3 while tp which determines 
a phase in the unitary operator remains an arbitrary pa- 
rameter. 
Thus the initial state of S can be determined as 

r x = V3(s z a z ), 

r y = V3(cos2tp(s z ) — sin2tp(a z )) , 

r z = V3(shi2ifj(s z ) + cos2ip(a z )) . (57) 

If we choose the phase tp = ir/4, we get 

r x = V3(s z a z ), 
r v = -v3(fr z ), 

r z = \/3(s z ). (58) 

Following the same line of arguments analogous to the 
random initial state of the assistant, a proper Hamilto- 
nian corresponding to the optimal determinant of 1/I2y3 

is 

H = —=a x s x + -(s y sma + s z ), 



(59) 



This Hamiltonian represents an Ising interaction between 
A and S together with an action of a magnetic field on 
the assistant. 



III. LIGHT AS AN ASSISTANT 

In this section we introduce another type of assistant. 
Specifically we show that the unknown density matrix 
of an ensemble of two-level systems (atom or spin) can 
be determined via interaction with a single mode of the 
electromagnetic field. In this case the type of the interac- 
tion Hamiltonian is fixed from the beginning. The atom- 
field interaction is studied within the Jaynes-Cummings 
model (JCM) 11 ; 12 . Thus the unitary time-evolution op- 
erator is known. We choose two commuting observables 
of the overall system and show that the initial state of 
S is linearly dependent on the expectation values of the 
two observables and their correlator. In this case, the un- 
known state of the spin can be characterized by repeated 
measurement of two commuting observables: the popula- 
tion difference of the atoms a z , and the photon number of 
the field a) a. This measurement supplies three averages: 
(<7 Z ), (a'a), and (a z a^a), which will be linearly related 
to the elements of the initial density matrix of the en- 
semble of the two-level atoms. (Note that since a z and 




a^a commute, (a z a^a) is recovered from the measure- 
ment data of a z and cva via the number of coincidences. 
This is similar to previous section, where (s z a z ) could 
be determined from measurements of s z and a z .) 



A. Interaction Hamiltonian 

The Jaynes-Cummings Hamiltonian for the interaction 
of a two-level system with a single mode field reads^ 



H = hujG z + hva}a + hg(a+a + <j_a') 



(60) 



where we restricted ourselves in the dipole approxima- 
te and v are the atom and field frequencies, re- 



tio: 
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spectively. a and a^ are field annihilation and creation 
operators, while &+ and &- are the raising and lowering 
spin operators, and g is the coupling constant. We shall 
denote 

A de f 

A = U! — V, 

for the detuning parameter. For our future purposes we 
note that A is a tunable parameter. 

The time evolution operator of the Jaynes-Cummings 
(JC) Hamiltonian, U(t), can be calculated in an exact 
manner [see— and references therein.]. 

In the eigenbasis of the two-level system U(i) reads: 

U(t) = e-^i^a+h) 

x cos[<v/^ + 5 2 ] - 1- l yJ-J? |+X+| 

V l V'P + r j 



_ ige -^(at a+ i)Bin[ty^T^] a|+)( _ | 

V + 9 2 
-ige-™ (* t »-i)™*^ a t|_ >(+ | 

+ e--^-)(cos^H-4^)|->(-|, 

(61) 

where |±) are the eigenstates of a z with eigenenergies 
E±. 

The unitarity of U(t) is satisfied because of the identi- 
ties 



t y/tf> + ,9 s 



\j$ + .9 2 



sin [t y/ip] 



r 



a cos 



tyjlf) 



(62) 



Having U(t) at hand, we can calculate any property 
of S + A, the overall system, we wish. We consider the 
most general form of the initial state for the atom. This is 
described by some general mixed density matrix p given 
by©. 



For the assistant, we shall assume that the single cavity 
mode starts its evolution from a coherent state with a 
known parameter a: 



M 2 /2 



E 

n=0 



1)1 



(63) 



where \a) is the eigenvector of the annihilation operator 
a, 

a\a) = a|a), 

and where \n) is the eigenvector of the photon number 
operator aJa, 



tion difference (a z )t reads 



2 oo 



(&z)t = yX!(' l + !)(c«+i - c„ 



n=0 



sin 2 (Q,„t/2) 
(^„i/2) 2 



■H^r.f; c„^^S- 



n=0 

oo 



+4 5 r y ^, 



71=0 



sin (0„t/2) 



^{Xn(t)} 



1- g 2 ^2(n+ l)(c„+i + c„) 



n=0 



sin 2 (Q n */2) 

(a.72) 2 



(69) 



(r a|n) = n|n). 

The assumption ([6"3")) on the initial state of the field is 
natural since these are the kinds of fields produced by 
classical currents^, and also, to a good approximation, 
by sufficiently intense laser fields. 

We assume the system and the assistant are initially sep- 
arated and do not interact with each other. As a result, 
the overall initial density matrix is factorized, 



V(0) = p\a)(a\, 



(64) 



where the initial state of the system in the eigen-basis of 
a z reads 



where r x ,r y , and r z are the unknown elements of the 
initial atom density matrix, which we want to find out, 
!K and 9 stand for the real and the imaginary parts of the 
argument in front of them, respectively. The parameters 
Xn{t), c n are defined as 



u\ de f 
X-n.it) = a 



n n t\ sin(ru/2) 

cos + iA 



i i n 



and 



d lf „~H 2 a2n 



(70) 



(71) 



\\+rA\+){+\ + {r x -ir y )\+){-\ 
+ (r x +ir tf )|-><+| + Q-r,)|->H. (65) 
The initial state of the field is given by 

oo oo n *m 

\a)(a\ = e-H 2 £ £ iL?Li n)(m |. (66) 
— n — n Vntyml 



n— TTi—0 



The state of the system + assistant at later time t can 
be calculated with the help of the unitary operator ([6l)l 



i>(t) = u(t)p(o)u\t). 



(67) 



Then the expectation value of any observable O of the 
overall system at time t is 



(6) = tr V{t)6 



(68) 



The next step is to calculate the two commuting observ- 
ables with which we can build up the initial state of the 
atom. Using (J5TJ) . (f6"5)) . ([6"6")l . and (|57|) . the atom popula- 



where the corresponding Rabi frequency, f2 n , is defined 
as 



def 



Q n = J V4(n+1).9 2 + A2. 



(72) 



The average number of photons in the cavity, (dta) t , can 
be calculated in a similar way 



(tfa) t = y^n-c„ 



n=0 



tE(' 1 + 1 )( c »+ 1 



-4gr x y^c» 



sm(n n t/2) ( 

~ - n 



sm 2 (n n t/2) 

(n n /2y 

HXn(t)} 



-4 S .E^^»- 



-'■■/; 



+g 2 r z yj(n + l)(c ra+ i +Cn) 



n=0 



sm 2 (n n t/2) 

(^n/2) 2 ' 



(73) 



The correlator of the two observables, (ova^a)*, which 
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FIG. 1: The dynamics of the determinant of the matrix M in the Jaynes-Cummings model for different values of 

the mean photon number in the cavity, n = 2, 5, 10 with two detuning parameters: A = WKHz and A = WOKHz. The figures 

are cited from reference-. 



amounts to the number of coincidences, reads 
(a z a*a) t = 



9 — ^(7i + 1) [(2n + 3)c„+i - (2n + l)cj 



Expectedly, these three quantities, i.e., the atom pop- 
ulation difference {cr z )t, the average number of pho- 
tons (fit a) 4, and the correlator of these two observables 
(<r 2 cra)t are linearly related to the three unknown pa- 
rameters r x , r y , r z of the initial atom density matrix: 



sin 2 (Q„i/2) 

' («„/2) 2 

-2gr x f>(2n+l) S ™ ( ""* /2) 3{ X n(*)} 



n=0 



i L n 



+2gr y JTc n (2n+l) S ' m{ ^ t/2) $t 



n=0 

(n + l)q 2 
r, > -I nc.„ x 



i Lin 



ra=0 



x [(2ra + 3)c n+ i + (2n + l)c„ 



sin 2 (fi„t/2) 

(n n /2r 



}■ 




B 



(75) 



The elements of the 3x3 matrix M. and the vector B 
are read off from Eqs. ([55]) - ([73)1 . They depend on 
the parameter a of the initial assistant state, on the de- 
tuning parameters A, coupling g of the JC Hamiltonian, 



(74) 



10 



and on the interaction time t. Thus, if the matrix Ai 
is non-singular, i.e., the determinant of M. is not zero, 
one can invert M. and express the unknown parameters 
of the initial atom density matrix via known quantities. 
Although the elements of M. are complicated, the deter- 
minant itself is much simpler. It takes the explicit form 

D{t) d = f dct[M] = 4A<? 2 e- 2 l Q l 2 

- ~ | a |2(n+m+l) 

x > > : — : (n — m) x 



n— m— 



n\m\ 
2 fn j- /o\ „:„ o + „:„2 



sin (Q n t/2) sin Q m t sin (fl m t/2) sinH, n t 






Q 2 fl 



(76) 



At the initial time t = 0, the determinant -D(O) is nat- 
urally zero, since the initial state of the overall system 
is factorized. Note that D(t) = for A = 0. Aquino 
et alr^ showed in the resonance case this scheme can be 
implied for the spin-^ systems but not for a quasi-spin-i 
system. Thus some non-zero detuning is crucial for the 
present scheme of the state determination of any two- 
level system. The value of D(t) changes by varying the 
detuning parameter A. It is seen in Figs. 1(a) 1(f) that 
for a non-zero detuning, A /0, the determinant (t) is 
non-zero for a certain period t > 0. (Obviously D(t) is 
zero when there is no photon in the cavity. ) On the other 
hand, large D(t) suggests that the state of the atom and 
the field are entangled^. 



Comparing figures Fig. 1(a) and Fig. 1(c) we see that 
although higher initial photon numbers n lead to bigger 
values for the determinant, they cause rapid oscillations 
in the value of the determinant. This makes the measure- 
ment process more difficult. (Note in this context that 
the determinant depends on the absolute value of a and 
n = \a\ 2 is the average number of photons.) 
If the average number of photons n = \a\ 2 in the ini- 
tial state of the field is sufficiently large, the determinant 
is nearly zero for intermediate times; see Figs. 1(e) and 
1(f) The reason for this collapse is apparent from (|T6")) 



and has the same origin as the collapse of the atomic 
population difference well known for the JCMi£. Each 
term in the right hand side of (|76[) oscillates with a dif- 
ferent frequency. With time these oscillations get out of 
phase and D(t) vanishes (collapses). However, since the 
number of relevant oscillations in D(t) is finite, they par- 
tially get in phase for later times producing the revival 



of D(t), as seen in the Figs. 1(e) and 1(f) 



It is seen that D(t) does not depend on separate frequen- 
cies bj and v of the two-level system and the field, only 
their difference A = uj — v is relevant. This is due to the 
choice of the measurement basis — see the left hand side 
of (j75|) — that involves quantities which are constants of 
motion for g — > 0. Comparing the figures Fig. 

and Fig. 



Fig. 1(b) Fig. 1(c) with Fig. 1(d) 



l(a)| with 
e)| with 



Fig. 1(f) one observes that the value of the highest peak 



of D increases by an order of magnitude when the detun- 
ing parameter changes from 10kHz to 100kHz. Note that 



in Eq. ([76]) for the determinant D(t) the contribution 
from the diagonal n = m matrix elements of the assis- 
tant initial state |a)(a| cancels out. Thus, it is important 
to have an initial state of the assistant with non-zero di- 
agonal elements in the {\n)} basis. 

Since the determinant D(t) is not zero for a realistic range 
of the parameters, the initial unknown state of the two 
level system can be determined by specifying the average 
atom population difference {(J z )t, the average number of 
photons (a*a)t, and their correlator (a z a)a)t- 



B. Random interaction time 

We saw in the previous sections that the success of 
the presented scheme is to a large extent determined by 
the ability to select properly the interaction time t, since 
this ultimately should ensure a non-zero (and sufficiently 
large) determinant D(t) since a small determinant will 
amplify numerical errors. 

We notice that the expectation value of an observable 
O, as it is described in section I, is the ensemble aver- 
age of an observable. This is mathematically denoted by 

tr pO . Now we have to take into account that the re- 
peated measurement of counting the number of the pho- 
ton in the cavity and the population difference of atoms 
are performed at a random t in each set of measurement, 
which obeys the Gaussian distribution. Thus we have to 
perform a time-average in the relevant time window as 
well. Since we just want to get a rough estimation about 
the consequence of such way of measurement on the value 
of the determinant, we avoid the tedious time- averaging 
calculation of (d^), (<r z ), and their correlation by making 
a shortcut and perform the time-averaging of the deter- 
minant itself. 

To quantify the robustness of the presented scheme it is 
reasonable to assume that there is no perfect control in 
choosing the interaction time. To this end let us assume 
that the interaction time t is a random, Gaussian dis- 
tributed quantity centered at to with a dispersion a and 
that an ensemble of measurements is performed to map 
out this spread. The corresponding probability distribu- 
tion P(t) of thus reads 



P(t) = — — e -(*-*°) 2 /(2<^ 

2na 



(77) 



Since n order to get a roughly estimation about the in- 
fluence of the presence of such randomness during the 
measurement procedure on the value of the determinant 
we make a short cut and time average the determinant 
D(t) over the same probability distribution 



< a t2(n+m+l) 



i\m\ 



D(t ) = AAg 2 e-^ 2 J2E- 

n— m— 

x (n - m) [w(Q, n , Q m ; t ) - w(tt m ,tt n ; t )} , 



(78) 
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(a) A = 10kHz, a = 0.1/j.s 



(b) A = 100kHz,a = O.Ol^s 



FIG. 2: The time averaged determinant D in the Jaynes-Cummings model as a function of to when the mean photon number 
in the cavity is n = 2, and g = 50kHz for different values of A and a; see Eqs. (|77|l - (|79|l . The figures are cited from reference^. 



where 



w(n n ,Q m ;t Q ) 



40 2 



{2e~ 



sinfto^r 



e 2 



e 2 



(n m +n n f sin[io ( fim + Q n) ] 



sin[(Q m - il„)t ]}- 



(O m -O n ) ! 



(79) 



It is seen that the oscillations of D(t) turn after averaging 

into exponential factors e _ ^"'° m±n "' ) and e~ £ 5~ n ™ in 
([751 IT9")) . due to which the averaged determinant D(to) 
gets suppressed for a sufficiently large "indeterminacy" er. 



This suppression is illustrated in Fig. 2(a) and Fig. 2(b) 
By comparing Fig. 2(a) and Fig. 2(b) we realize that 



when the dispersion a grows by one order of magnitude, 
the value of the averaged determinant drops dramatically. 



IV. CONCLUSION 

It is important to implement the single-apparatus to- 
mography for a situation with a physically transparent 
measurement base and with a realistic system-assistant 
interaction. Here we carried out this program for a two- 
level atom (system). The non-commutating elements of 
the density matrix of a quantum system can be deter- 
mined by simultaneous measurements of commuting ob- 
servables. This can be done by introducing another sys- 



tem which its initial state is known, called assistant. The 
two systems interact and then by performing repeated 
measurements of one observable belonging to the system 
of interest and the other one to the assistant, one is able 
to count the events in repeated experiments. This yields 
to a one-to-one correspondence between the initial den- 
sity matrix and the collected data. We displayed several 
examples. In the first scries of examples we considered 
the assistant to be a two-level system with a known ini- 
tial state. We showed that the full initial density matrix 
of the system can be determined via simultaneous mea- 
surements of the occupation probabilities for the energy- 
levels of two systems. We discussed two different initial 
states for the assistant and showed that this procedure is 
feasible even when the assistant starts its evolution from 
completely disordered state. 

The dynamical processes which properly yield the deter- 
mination of the initial state of the system can be de- 
scribed by an Ising or Heisenberg interaction Hamilto- 
nian, regarding to the situation. 

We also showed that another way of reconstructing the 
state of a two-level system is to let it interact with a 
specific type of environment, namely, a coherent single 
mode of an electromagnetic field. In this case the interac- 
tion Hamiltonian is described by the Jaynes-Cummings 
Hamiltonian. In this sitation, the measurement of the 
simplest set of observables related to the energies of the 
atom and field yield the determination of the unknown 
initial state of the atom. 
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